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Abstract. We consider algebraic functions z satisfying equations of the form 
(0.1) ao2" + aiz"'^ + a22"^ + • • • + a«^™" + a«+i = 0. 

Here m > m\ > . . . > to„ > 0, to, to^ e N, and z ~ z{ao, . . . , a„+i) is a function 
of the complex variables oq, . . . , a„+i. Solutions to such equations are classically 
known to satisfy holonomic systems of linear partial differential equations with 
polynomial coefficients. In this paper we investigate one of such systems of dif- 
ferential equations which was introduced by Mellin. We compute the holonomic 
rank of the Mellin system as well as the dimension of the space of its algebraic 
solutions. Moreover, we construct explicit bases of solutions in terms of the roots 
of (|U.1|) and their logarithms. We show that the monodromy of the Mellin system 
is always reducible and give some factorization results in the univariate case. 



1. Introduction 

Consider the class of algebraic functions which satisfy algebraic equations with sym- 
bolic coefficients, i.e., equations of the form 

(1.1) aoz"" + aiz""' + a2z"'^ + ... + a^z'^- + a„+i = 0. 

Here m > mi > . . . > m„ > 0, m, G N, 2; = z{aQ, . . . , a„_|_i) is a function of the 
complex variables oq, . . . , ctn+i- Any solution to (jl.lj) satisfies the A-hypergeometric 
system defined by Gelfand, Kapranov and Zelevinsky jS] (see also OIII]) associated 
to the matrix 

'1 1 ... 11 

m nil ■ ■ ■ f^n 

and the homogeneity vector (0,-1). In particular, a solution to (jl.lj) has a double 
homogeneity property and can therefore be considered as a function of n variables. 
This implies that we can arbitrarily prescribe the values of any two nonzero coef- 
ficients in fll.lj) without losing any essential information on the general solution to 
this equation. It turns out to be convenient to divide p.ip by — a„+i and then set 
y = {—ao/an+iY^^^z, reducing p.lj) to an equation of the form 

(1.2) + xiy'"^ + . . . + x^y™" - 1 = 0. 

A classical result of Mellin from 1921 (see JUl) states that solutions to such equations 
satisfy certain systems of partial differential equations of hypergeometric type. More 
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precisely, if y{x) = y{xi, . . . , x„) is a solution to ()1.2j) . then it satisfies the following 
system of n partial differential equations: 



m,-l m-l 



Y\ (miOi + . . . + rUnOn + mk + 1) JJ^ {m[9i + . . . + m'^On + mk - l)y{x) = 

k=0 k=0 

(1.3) i-^r^rn-^^, j = l,...,n, 

where 9j = Xj-^ and m- = m — mj. In this article we will only consider algebraic 
equations of the form (|1.2p . The system (jl.3p will be referred to as the Mellin system 
of equations associated with il.^) . 

The goal of this paper is to describe all the solutions of the Mellin system of 
equations and to construct explicit bases of local complex analytic solutions for any n 
and any value of the greatest common divisor d of m, mi, . . . , m„. We describe the 
space of algebraic solutions and we exhibit explicit non-algebraic solutions in terms 
of the roots of and their logarithms. These problems arose in the framework 
of the general theory of hypergeometric functions and their relations with algebraic 
equations, see |3j , [El , [SI , (JHI; and [T7] . 

The special case of the trinomial equation y"^ + xy^ — 1 = (which corresponds 
to ordinary Mellin equation) was thoroughly studied in The Mellin equation is 
always reducible and unless p = m — 1, it is never the linear differential equation 
with polynomial coefficients of smallest order satisfied by the roots of the trinomial 
equation. Moreover, in this case all solutions to the Mellin equation are algebraic 
except in the case when d = 1 and p < m — 1, where we exhibit an explicit non- 
algebraic solution, as a special case of our main Theorem 14.31 

We show that the monodromy of the Mellin equations is always reducible. This 
could be deduced when d = n = 1 from the general results in [T7j since in this 
case the homogeneity (3 has integer coordinates. In the last section, we give some 
factorization results in the univariate case. 

Acknowledgements: The authors are grateful to Martin Mereb for useful com- 
ments on Theorem 14.31 and to Michael Singer for interesting discussions and for his 
help in the proof of Proposition 15.51 T. Sadykov is greatly indebted to his col- 
leagues and staff at the Max-Planck Institut fiir Mathematik in Bonn for creating a 
wonderful working atmosphere. 

2. Convenient bases and generating solutions of the Mellin system 

Mellin calls the solution of ()1.2j) which assumes the value 1 for xi = . . . = x„ = 
the principal solution of this equation. He shows (see [TU]) that the power series 
expansion of the principal solution to (jl.2j) around the origin is given by 



(2.1) y,r{x)= 



Here |z/| = z^i + . . . + z/„ and the empty product is defined to be 1. 
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For a multi- index / = (ii, . . . , and a complex vector x = (xi, . . . , x„) G C" we 
use the standard notation = x^^ . . . x^". By B we denote the set 

B = {(ii, . . . , i„) G Z" : < < m — 1, for any j = 1, . . . , n}. 

Theorem 2.1. T/ie dimension of the space of analytic solutions to at a generic 
point, i.e., the holonomic rank of the Mellin system, equals rrf'. There exists a 
basis {//}/G-B around the origin such that 

(2.2) /,(x)=a;7/«,..-,0, 
where fi is analytic at and //(O) 7^ for any I E B. 

We say that a system of functions of the form ()2.2j) forms a convenient basis in 
the solution space of the Mellin system ()1.3|) . 

The proof of this theorem can be deduced from Theorems 2.8 and 3.1 in ^1] by 
adjusting the arguments in [H] to the case of a nonconfluent hypergeometric system. 
For the benefit of the reader we give now a proof which uses the results of P3| as 
little as possible. We will reprove a part of this theorem in the next section, as 
a consequence of the relation with the A-hypergeometric system associated to the 
exponents m, mi, . . . , m^, 0, and its translation to a Horn system |P. 

Proof. For the sake of brevity we will use the notation M = (mi, . . . ,mn), M' = 
(m^, . . . , m„), 9 = {9i, . . . ,9n) and denote by (■, ■) the scalar product. We multiply 
the j-th equation in p.3|) with x™. This will not affect the space of holomorphic 

solutions to p.3|) . Let us denote by Pj{9) the polynomial YTk=o^i{^^^) + + 
1) Y[k=o ((^ ' ^) + ~ by Gj the differential operator 

(2.3) G, = xYP,i9) - (-l)-.m-xf j = 1, . . . ,n. 

Let V be the Weyl algebra in n variables 0, and define A4 = V/ X]j=i ^e 
the left P-module associated with the system of differential operators ()2.3p . Let 
R = C[zi, . . . ,Zn] and R[x] = R[xi, . . . , x„] = C[xi, . . . , x„, zi, . . . ,Zn]. We make R[x] 
into a left P-module by defining the action of dj on R[x] by 

(2.4) = 

Let us define the operators Dj : R[x] R[x] by setting 

d 

(2.5) Dj- = 2;^ — + XjZj, J = l,...,n. 

It was pointed out in pP that the operators ()2.5|) form a commutative family of 
D-linear operators. Let D denote the vector {Di, . . . , Dn). For any j = 1, . . . , n we 
define the operator Vj : R[x] — > R[x] by Vj = zJ^Dj. This operator commutes with 
the operators dk since both Dj and the multiplication by zj'^ commute with dk- 
Moreover, the operator Vj commutes with Vfc for all 1 < j,k < n and with Dk 
for j 7^ /c. In the case j = k we have VjDj = Vj + DjVj. 
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Let US define a family of linear operators {Wj}""^-^ acting on the ring R[x] by 



m— 1 



k=0 

We claim that the family of operators {Wj}^^^ is commutative. Indeed, using the 
fact that {Dj}'j^i is a commutative family of operators and the identity V™-Dj = 
{Dj + m) V™" we conclude, that 

W,Wk = P,{D)Pk{D,, ...,D,+m,...,Dn) VfV^ = 
Yl ((M,D)+m/i + l) Yl ((M',D) +m/i- 1) 

y At=0 /x=0 y 

Computing the composition of the operators Wj and Wk in the reversed order, we 
arrive at the same expression. 

Using Lemma 2.2 and Theorem 2.4 in ^3] we conclude that the P-module asso- 
ciated with the modified Mellin system ()2.3p is isomorphic to the quotient 

R[x] j W,R\, 

By the commutativity of the operators Wj and Lemma 2.7 in the dimension 
of this quotient as a complex vector space equals the product of the degrees of the 
operators in the Mellin system, i.e., m". This justifies the first claim of the theorem. 

Let us now prove the second statement of the theorem. Let / G -B and consider 
a power series of the form 

(2.6) /,(x)= ^(^)^'' 

where mNp + / = {(si, . . . , s„) G Nq : = mp^ + i^, for some G Nq}. To compute 
the action of the operator Gj on //(x) we notice that 

(2-7) -r^ = n(^^-^) 

and that 

Since nfc!ro^(^i — ^) = for any I E B and any j = l,...,n, it follows, that 
Gjfi{x) = if and only if 

m— 1 

(2.8) <^is)Pj{s + !) = (-l)™^m"'(^(s + mcj) JJ {sj + m + ij-k). 

k=0 

Here ej = (0, . . . , 1, . . . , 0) with 1 in the j-th place. It is always possible to satisfy this 
relation by choosing {p{s) in an appropriate way. Indeed, the product nfc!ro^('^i + 
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m + ij — k) is never zero and hence the recurrent relations ()2.8|1 can always be 
solved. Thus for any I E B there exists a formal series solution to (jl.3|) of the 
form ()2.6|) . Notice that the resultant of the principal symbols of the operators in 
the Mellin system ()1.3p is different from zero at the origin. By Proposition 2 in ^2] , 
any formal series solution to the Mellin system of the form ()2.6|) converges in some 
neighbourhood of the origin. The m" series {fi{x)}i(zB are linearly independent 
since their initial monomials are different. Thus they form a convenient basis in 
the space of analytic solutions to the Mellin system. This completes the proof of 
Theorem 12.11 □ 

Definition 2.2. We call an analytic solution y{x) to the Mellin system ()1.3|) defined 
near the origin a generating solution to this system if for any I E B, 

djy{0) := . ^ 0. 

dx'^ . . . 



As we will see in Theorem 12. 4[ a generating solution gives rise to a basis in the 
solution space of the Mellin system. On the other hand, the sum of the elements of 
a convenient basis is a generating solution to p.3p . 

Example 2.3. The Mellin equation associated with the quadratic equation y"^ + 
xy — 1 = has the form 

(2.9) (x^ + 4)y" + xy' -y = 0. 

A convenient basis in the solution space of this second-order differential equation 
is given by the functions = x, f2{x) = + 4. Another basis is given by 

the roots yi = |(— /i ± f2) of the quadratic equation y"^ + xy — 1 = 0, which are 
generating solutions. 

We denote by e the m-th primitive root of unity e := e^'^*/™. 

Theorem 2.4. Let y{x) be an analytic function near the origin which is a solution 
to the Mellin system M.'J\) . Then, y{x) is a generating solution if and only if the 
family of functions yi{x) = y{e^'^Xi, . . . , e^"Xn), I = {ii, . . . , in) ^ B is a basis in 
the solution space of / ti.,'y)) . 

Proof. Assume that y{x) is a generating solution. The differential operators 6^ = 
^'^af^ 'B^ invariant under the change of variables tk = s'^^Xk for any k = 
1, . . . ,n. Indeed, Xk-£^ = ^~**^A:^^ = '^^^ invariance of follows from 

the identity = 1. Since the differential equations in the system ()1.3|) only contain 
such operators, it follows that the function yi{x) = . . . solves ()1.3|) 

for any I E B. 

Let us show that the family of m" functions {yi{x)}i(zB is linearly independent. 
This will prove the assertion of the theorem. Let 

y(^) = 5Z v{s)x' 
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be the Taylor series expansion of y{x) around the origin. For a given I & B consider 
its subseries defined by 

(2.10) fj{x)= Yl ^(^)^^ 

where mNg + / = {(si, . . . , s„) G Ng : Sfc = mpk + ik, for some pk G No}. By assump- 
tion y{x) is a generating solution to Hence the family of series {fi{x)}iizB is 
linearly independent since their initial monomials are different. Being subseries of a 
convergent power series, they all converge in some neighbourhood of the origin. It 
follows directly from the definition of fi{x) that y{x) = J2i£B fii-'^)- 

Suppose that the family of functions {yi{x)}i(zB satisfies a linear relation 

(2.11) J2ciyi{x)=o 

for some c/ G C. Let J = (ji, . . . , j„) G -B be a multi-index. Using the definition of 
the function fj and the identity = 1 we conclude that 

fj{e''xi,...,e'-Xn) = Yl ^is){e''xiy' ...{e'-XnY" = 

s&mNo+J 

J2 e''^'+-^'"'"ip{s)x' = £^^'^V(s)^' = ^^^'^Vj(^), 

where (/, J) denotes the scalar product. Thus the relation 1)2.111) can be written in 
the form 

= ^ ciyi{x) = Y civie'^xi, . . . , e'"Xn) = 
leB i£B 



(2.12) 5^c,5^/j(e^^a:i,..., 

/GB JG-B 




The assumption that y{x) is a generating solution to the Mellin system ()1.3|) implies 
that the family of functions {/j}jg_b is linearly independent. It therefore follows 
from (E:T^ that 

(2.13) Y c/^^'^"^^ = 0, for any J e B. 

I<^B 

The matrix of this system of linear homogeneous algebraic equations is the n-th 
exterior power of the Vandermonde matrix 





1 


1 


1 


1 


\ 




1 


e 


. 








1 




. 


^2(m-l) 




v 


1 




^2(m-l) 


^{m-l)2 


/ 



This Vandermonde matrix is nondegenerate since the numbers 1, e, e^, . . . , ^ are 
all different. It is known that the eigenvalues of the exterior product of two matrices 
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are products of the eigenvalues of the factors. Thus the exterior product of two non- 
degenerate matrices is nondegenerate and hence so is the matrix of the system (j2.13|) . 
This imphes that Cj = for any I E B and hence the functions {yi{x)}i(zB are hn- 
early independent. 

Assume now that y{x) is not generating and let J G -B such that djy{0) = 0. 
Then, djyi{0) = for all I E B. By Theorem 12.11 there exists a solution fj with 
djfj{0) = 1. Then, the family {yi{x), I G B} cannot be a basis of solutions of the 
Mellin system. This completes the proof. □ 

Recall that we denote d = GCD(m, mi, m2, . . . , 

Corollary 2.5. If d > 1 then {ypr{e^^Xi, . . . ,e''"Xn)}i<^B is a basis in the solution 
space to the Mellin system / li.,^) . Here ypr{x) is the principal solution h2.1\] to the 
Mellin system. In particular, under the above assumption all solutions to the Mellin 
system are algebraic functions. 

Proof. \{ d> 1 then the term miVi + . . . + ninVn ~ + 1 in the numerator of the 
coefficient of (j^.lj) is never zero. In particular, it is nonzero for [ui, . . . , G B and 
hence ypr{x) is a generating solution to the Mellin system. The conclusion of the 
corollary follows now from Theorem 12.41 □ 

Example 2.6. Consider the algebraic equation 



A generating solution for the corresponding Mellin system is given by the function 



where u{xi, X2) = y 108 + 86x1X2 — 8xf + 12^^81 + 54xiX2 — 12xf + 12x1 — Sxfx^. 
A simple computation in Maple shows that all of the initial exponents for ()2.14|) are 
present in the series expansion of ?/(xi, X2), as predicted by Corollarv l2.5l In contrast 
with Example 14.21 below, in this case it is possible to give a generating solution to 
the Mellin system associated with ()2.14j) in terms of the solutions to ()2.14|) itself. 

Remark 2.7. From the point of view of the original algebraic equation ()1.2|1 it is, 
of course, very unnatural to make the assumption that d > 1. Indeed, the change of 
unknown u{x) = y{xY reduces ()1.2|) to an equation with relatively prime exponents. 
However, it turns out that the Mellin system associated with this transformed equa- 
tion differs essentially from the Mellin system associated with the original algebraic 
equation. This is illustrated by Examples 15 . 71 and 15 . 3[ and by Examples 12 . 61 and 14. 2[ 
since in case d > 1 all solutions are algebraic by Corollary 12. 5^ while as we will see, 
in case d = 1, n = 1, mi <m — lor(i=l,n>l, there always exists a non- algebraic 
solution. Thus from the point of view of the Mellin system the case when d > 1 
cannot be trivially reduced to the case of relatively prime exponents. 



(2.14) 



/ + Xiy^ + X22/^ -1 = 0. 




(3.1) 



A 
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3. A-HYPERGEOMETRIC, HORN HYPERGEOMETRIC AND MELLIN SYSTEMS 

Given the equation its roots (thought of as functions of a = ao, . . . , fln+i)) 

satisfy the ^-hypergeo metric system with parameter (3 := (0,-1) (cf ^IGJ), where 

1 1 ... 1 1 

m nil ■ ■ ■ iTT'n 

Namely, it is the left ideal in the Weyl algebra C[ao, . . . , a^+i, do, ... , dn+i] generated 
by 

(3.2) the toric operators 9" — for u,v G N"^^ with A ■ u = A ■ v, 

n+l n 

(3.3) and the Euler operators ''^^^Ljdj and moo + "^^rrijajdj + 1. 

j=0 j=l 

Note that the corresponding A-hypergeometric system is the same as the hyper- 
geometric system associated with the matrix 

1 1 ... 1 1 

m/d mi/d ... rrin/d 

and homogeneity /?' = (0, —1/d). 

Consider now the following matrix B E 

/ —nil —^2 ■ ■ ■ \ 
m ... 
m ... 



(3.4) 



A' 



(3.5) 



B 



... m 

\ -m; -m'2 . . . -m; / 

where as above m[ = m — nii, i = 1, . . . , n. Let c = (— 1/m, . . . , 0, 1/m) and denote 
by HornB(c) = (Hi, . . . ,Hn) the corresponding Horn system, as in Definition 2.1 
in jH], which we now recall. For j = 1, . . . ,n, the Horn operator Hj in the variables 
w = {wi, . . . , Wn) equals 
(3.6) 

m-l ruj-l ^ m.-l ^ 

Y\ {rnej-k)~Wj TT (-mi6'i-. . .-m„6'„, k) TT (-m'-L6'i-. . .-m|^6'„H k). 

k=0 k=0 k=0 

We deduce from Lemma 5.1 and Corollary 5.2 in j6j the following result. 

Lemma 3.1. Let '0(^1, . . . , w„) be a function ofn complex variables and let ip{a) be 
the following function of n + 2 variables ( defined in suitable simply connected open 
sets): 



(3.7) 



mi "'1 
"0 "n+l 



"0 "n+l , 



Then, (p is a solution to the A-hypergeometric system with parameter (3 if and only 
if if) is a solution to Horng(c). 
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Setting ao = 1, a„+i = —1, we deduce the relation 

^(1, a,,..., -1) = (-1)^ z^((-ir'iar, . . . , (-ir^a^T)- 

We thus make the change of variables Wj = (— and we get that given a 
holomorphic function ipi^Wi, . . . ,Wn) as above, the function 

(3.8) /i(xi, . . . , x„) := ^{w{x)) = ^{{-ir'^xT, • • • , (-1)"^"0 
verifies 

(3.9) (-l)™yu(ai, . . . , a„) = ^^(l, ai, . . . , a„, -1). 
We now translate the Horn operators in the w variables to the x variables. Since 



^isS- = ^^ia!-' get for all J = 1, ... , m, 

m-l rrij-l . 

H- ■■= \{{e, -k)- i-ir'^x- n ( --^1 0n-- 
J ±± J J ±± \ mm 

k=0 k=0 ^ 



„ m„ „ 1 



(3.10) n 

\ m mm 
k=0 ^ 

where now 6'i, . . . , 6'„ refer to derivatives with respect to the variables (xi, . . . , Xn), 
i.e., 6i Xi i 1, . . . , 72.. 

Note that 11^=0 (^i -k) = xf (^j . Multiplying H'- by (-l)'"+im'^ and divid- 
ing by X™ we get the Mellin operators ()1.3p . 
We immediately deduce: 

Lemma 3.2. A holomorphic function fj.{xi, . . . , Xn) is annihilated by the Mellin op- 
erators il.!^) if and only if the function ijj{wi, . . . , Wn) verifies the Horn system ( I,?. 
where fi and ip are related by I'J. (in suitable open sets). 

We now take advantage of these translations to provide a new proof of the rank of 
the Mellin system stated in Theorem 12.11 Note that the greatest common divisor g 
of the maximal minors of B equals g = m'^~^d. The normalized volume of A 
equals m/d. Moreover, by the results in the lattice ideal Ib is a complete 
intersection. We deduce from jHl Theorem 2.5] that the rank of HornB(c) equals 
m^~^d X m/d = m". Using the translation of solutions in Lemma (3. 2t it is easy to 
check that the holonomic rank of the Mellin system is also m". 

Theorem 3.3. In suitable open sets, given an A-hypergeometric function <f{a) with 
parameter f3 = (0, —1), the function of n variables defined by specialization 

(3.11) /i(Xi, .. .,Xn) = <^{l, Xi,..., Xn, -1) 

is a solution to the Mellin system. 

Proof. Given a choice of m-roots of the coordinates, define 
ipiwi, ...,Wn):= l^ie ""'w^ ""-w;. 
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j_\m 



where e' is a primitive m-root of —1. Note that Wi = (— l)'"^ (^e' "^^wl 

1, . . . ,n. By Lemma I3.2[ it is enough to verify that ip is a solution to the Horn 
system. Now, by Lemma f3.H this happens if and only if the function (p' defined by: 



satisfies the A-hypergeometric system with parameter /3. But calling A an m-th root 
of — cio/'^n+i we deduce from the definitions of the functions and the homogeneity 
(0, —1) satisfied by (p that 

= ^-V((-^)A--,...,(-^)A--) 

= A-V(l, — ^A-^S . . . , -^A-™", -1) = ip(a), 

and we are done. □ 



4. Solutions in terms of roots 

Mellin not only observed in ^lOj that the roots yi{x), . . . ,ym{x) of the algebraic 
equation ()1.2|) 

I/™ + xiy""' + ... + x„?/™" -1 = 

satisfy the Mellin system ()1.3|1 . but he also made the following easy observation. 
Given the principal solution ypr{x) in ()2.1|) . all the solutions have the form 

riypriri"''xi,...,ri'^"Xn), 

where rj runs through the m-roots of 1. 

It is also clear that for any choice of / = (zi, . . . , z„) G N", the function yi{x) = 
ypr{s^^xi, . . . , £:*"x„) is a root of the algebraic equation 

(/) y'^+e'^xiy'''^+. . .+e'"Xny'^" -I = 0, 

where, as in Section|2l we denote e = e^'^*/"^. Moreover, as we have already remarked, 
all these functions lie in the solution space of p.Hj) . So, the Mellin system has not 
only the roots of the algebraic equation p.2j) as solutions, but also the roots of the 
associated equations (/). It is clear that we only get m" different equations this 
way, and we could parametrize them taking I E B. We denote by Y the C-vector 
space generated by all the roots of all the equations of the form (/). The roots of 
the original equation ()1.2j) correspond precisely to the functions y[jm-i,...jmn)) j — 
0, . . . ,m — 1. Note that these m functions are all distinct, but of course they are 
not linearly independent if c? > 1 or if = 1 and mi < m — 1. In case c? > 1 we 
have seen in Corollary 12.51 that Y coincides with the solutions space of the Mellin 
system. 



Example 4.1. The polyquadratic equation. Consider the equation 
(4.1) y'^'' + xy'' - 1 = 0. 
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Its 2k roots are y {—x ± V x'^ + 4)/2, j = 0, . . . , A; — 1, where e is the primitive k-th 
root of unity. Observe that these roots span a vector space of dimension 2. A basis 
in the 2fc-dimensional solution space of the Melhn equation associated with 1)4.11) is 

given by the functions ^/—x ± Va;^ + 4e^, j = 0, . . . , k — 1. 

For the rest of this section we will then assume that d = 1. Our task is to explain 
and generalize the following example. 

Example 4.2. The general cubic. Consider the algebraic equation 
(4.2) + xiy^ + - I = 0. 

The corresponding Mellin system is given by 



279i{9i-l){9i-2)y = xl{29i + 02 + l){29i + 9^ + 4){9i + 29^ - l)y, 
-2792(92 - 1) {92 - 2)y = xl{29i + 92 + l){9i + 292 - l){9i + 292 + 2)y. 



The three solutions to ()4.2|) are linear combinations of the following functions: 



yi{xi,X2) = xi, 1/2(2:1, X2) = ^/u{Xi,X2), y3iXi,X2) 



2 



3x2 — xf 
\/u{xi, 



X2 



where m(xi, 0:2) = 108 + 36x1X2 -8x^ + 12^/81 + 54xiX2 - 12x? + 12x1 - 3a;?x^. The 
expansions of y2 and y^ into power series around the origin are given by 

1/2 = 6 + -X1X2 — ^^1 + ~ ^^1^2 + • • • ; Z/3 = — -Xi — — a;iX2 + . . . 

In this example, the set of initial exponents is given by i? = G : < z, j < 

2}. The expansions of yi, y2, and 2/3 into power series contain 7 of the 9 elements in B. 
This could be also checked differentiating ()4.2|) implicitly and evaluating at 0. The 2 
missing exponents are (2, 1) and (0,2). Thus the solution space of ()4.3p contains a 
two-dimensional subspace which is not spanned by solutions to ()4.2|) . 

Remark that the roots of the equations (/) above with any index set of the form 
Ij = (jmi, . . . , jm„) for any j = 1, ... ,m — 1, are just constant multiples of the 
roots of the original equation p.2|) (corresponding to j = 0). More generally, there 
are ni^'^ subsets G^^\ . . . ^G^^" •* of algebraic equations of the form (J), each of 
them consisting of m equations, such that for any pair of equations (J), (J) in the 
same subset G^''^ , the corresponding solutions only differ by multiplicative constants. 
Denote by T = /('^""')}, where each I^^^ e B is a choice of an index set in 

the subset G^'^^ of algebraic equations, for all A; = 1, . . . ,m"~^. The analytic roots 
of the algebraic equation (J*^^-') defined in a small neighborhood of the origin will be 
denoted by yf\ ...,ym. 

Given a complex vector c = {c^^\ . . . , c*^™"" ^), denote by Xc the function defined 

by 

fe=i 
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where log denotes a holomorphic branch of the logarithm defined in the union of the 
images of the roots. 

When the exponents m, mi, . . . ,m„ are relatively prime, is it shown in |3J Th. 
2.4] (see also Cor. 3.6]) that the dimension of the vector space generated by the 
solutions to the algebraic equation ()1.2|) is equal to m if mi = m — 1 and it is equal 
to m — 1 if mi < m — 1. Thus, when d = 1 the only possible linear relation between 
the solutions to the algebraic equation ()1.2p is the obvious one: the roots of ()1.2|) 
sum up to zero if the coefficient by y"^"^ in (|1.2j) is zero. 

Our main result is the following. 

Theorem 4.3. The space Y coincides with the space of algebraic solutions to the 
Mellin system of equations around the origin. Its dimension equals dimY = m" — 
jjfi-i ^ I case mi = m — 1 and dimF = m" — m"~^ otherwise. The space of 
relations 

k=l 

has dimension dimi? = m"'~^ — 1 in case mi = m — 1 and dimi? = m"~^ otherwise. 
Moreover, for any nonzero c E R, the function Xc defined by ^.4\ ) is a non- algebraic 
solution to the Mellin system, the space 

S:={xc/ceR} 

has dimension dimS" = dim/2 and the sum S(BY is direct and equals the full solution 
space of the Mellin system. 

When n = 1 , S" is nonzero only when mi < m — 1 and in this case it is spanned by 

the non-algebraic function x = Vi log(l/i) H \-ym log(i/m), where ?/i(x), . . . , ym{x) 

denote the roots of y"^ + xy"^^ — 1 = 0, while the space of algebraic solutions Y has 
dimension m — 1. This is the content of Theorem 3.5 in |4j for A-hypergeometric 
systems associated to monomial curves, which can be directly translated to the 
Mellin equation via our computations in Section 01 So, Theorem 14.31 is inspired by 
and generalizes this result. When translating the Mellin system to the homogeneous 
setting as in Section El the ocurrence of logarithmic solutions is explained by the 
results in fS]. 

Example 4.4. (Example 14.21 continued^ In this case dimi? = dim 5* = 3 — 1 = 2, 

as we have seen. There are 3 subsets of 3 equations each and we could choose the 
following representatives to form our set F: 

(Ji) y^ + xiy"^ + X2y -1 =0, 

(/2) + xiy'^ + ex2y - 1 =0, 

ih) y^ + xiy'^ + e^X2y -1 = 0, 

where e is a primitive cubic root of 1, and Ji = (0, 0), /2 = (0, 1), /a = (0, 2). 

A basis of relations in R is given for instance by the two relations (1,-1,0), 
(1,0,-1), since the sum of the roots of each of the 3 equations above equals the same 

(k) 

value xi. On the other side, the space Y spanned by the functions yj , j, = 1, 2, 3 
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has dimension 7. So a basis of the solutions near the origin is given by choosing 
seven of these roots together with the two non-algebraic functions X(i,-i,o)! X(i,o,-i)- 

Proof of Theorem \4 'A We begin by computing the dimension of Y . Call (p^, for 
u G Ng, the coefficients in the expansion ()2.H) of the principal solution y = ypr of 
the algebraic equation ()1.2|) : 

(-1)'''' nLl"l^("^l^l + • • • + "^ri^^n - m/i + 1) 

= ^jTT^^j • 

Note that when d > 1, the solution is generating. This is proven in Theorem 12 .41 
whose proof we now generalize with similar notations. 

Let //,/ G B, be defined as in ()2.10p . Observe that when ipi = 0, then ipj' = 
for all /' G / + mNg. It follows that when ipi = then // = 0. Let S' = {/ G 
B : ifi ^ 0}. We claim that dimF = #5'. Clearly, y = J2ieB' fi family 
{// : / G B'} is linearly independent. Suppose that the family of functions {yi}ieB 
satisfies a linear relation 

(4.6) ^cm(a;) = 0, 

for some C/ G C. Then, 



leB 



JGB' \l£B J 

which implies that 

(4.7) J2 (^i^^^"^^ = 0' any ^ 

/GB 

In fact, {ci)i(zB satisfy ()4.6p if and only if these last equalities hold. Since the full 
matrix {e^^''^'^)jj^B is of maximal rank, it follows that the system of equations (|4.7p 
has rank ^B' and so the dimension of the space of linear relations R' among the 
algebraic functions {yi, I G B} is equal to #i? — #5'. Dualizing, we conclude that 
their linear span Y has dimension ^B', as wanted. 

We now prove that the cardinality of B' equals m" — m""^ + 1 in case mi = m — 1 
and dimF = m" — m"~^ otherwise. Call B" = B \ B' its complement in B. It 
is clear that u G B" if and only if there exists a positive integer yu strictly smaller 
than 1 1^1 such that 

miUi H h m„z/„ - m/i + 1 = 0, 

which we simply write as (M, u) — mfi + 1 = 0, where M = (mi, . . . , m„). Note that 
when this equality holds, it follows that the integer /i satisfies 1 < /i < since 
< rrin < ■ ■ ■ < mi < m and all Ui are non-negative. Moreover, it is easy to check 
that the equality /i = |z/| can only hold when mi = m — 1 and fi = ui = 1 and uj = 
for any j > 1. So it is enough to prove that the cardinality of the set 

{u e B : (M, u) = -1 mod (m)} 

equals m"~^. We claim that for any r = 0, . . . , m — 1 there exist precisely m""^ = 
m"/m solutions z/ G -B to the modular equation {M,i>) = r mod (m). Since we 
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are assuming that d = 1, there exists an integer n-tuple X = (xi, . . . , x„) satisfying 
(M, X) = 1 mod (m) and so the equation (M, v) = r mod (m) has a solution 
for any r. Then it is easy to see that for any z/, {(M, z/) mod (m), (M, z/ + X) 
mod (m), . . . , (M, z/ + (m — 1)-^) mod (tti)} takes all the values r = 0, . . . , m — 1 
precisely once. Moreover, two such subsets are disjoint or equal, from which the 
result follows since the cardinality of B is m". 

It follows that the space R' of linear relations among the algebraic functions 
{?//, / G -B} has dimension m"~^ — 1 in case mi = m — 1 and m*^"^ otherwise. The 
relations R define in principle a subspace of R' corresponding to those vectors c 
which have a special pattern of repeated coordinates, but it is easy to see that it 
has the same dimension. Indeed it is very easy to find a linear relation between the 
sum of the roots in each equation (which equals minus the coefficient of y™"^, which 
is a constant multiple of xi if mi = m — 1 and otherwise), with the sum of the 
roots of the original equation ()1.2|) . This gives m"~^ — 1 independent linear relations 
in R and we get moveover the relation c = (1, 0, . . . , 0) in case mi < m — 1. We 
can exphcitly hst a basis of the space R\ the vectors {ei — 62, ei — Cmn-i} when 
mi = m — 1 and the vectors {ci, em^-i} otherwise. 

Then, R = R' and we recover in particular the quoted result in |31 Th. 2.4] about 
the only possible linear relation among the roots of a single sparse algebraic equation 
with generic coefficients. 

In order to see that for any c G i? the function Xc is a solution to the Mellin system, 
we will need the translations described in Section El between the Mellin system and 
the corresponding Horn and A-hypergeometric systems. This will allow us to use 
the results in |1] based on the study of the relation between lattice hypergeometric 
ideals and A-hypergeometric ideals explained in jB] . We sketch the main ingredients 
of the proof and refer the reader to |31 IH] for the details. Since the lattice generated 
by the columns of the matrix B G z^"'~^^^^"^ in ()3.5|) equals m"~^ and the associated 
lattice ideal is a complete intersection, the Horn system ()3.6|) gets translated by 
homogenization (as in Lemma [3.11) to the hypergeometric system Tijs corresponding 
to this lattice. We get m"~^ hypergeometric systems Tij, i = 0, . . . , m"~^ — 1 which 
are torus translates of the A-hypergeometric system, all with the same homogeneity 
/3 = (0, —1) and such that all solutions of any Tii satisfy the system Tig. The different 
functions yj which satisfy the Mellin system, and thus the Horn system after the 
change of variables in Lemma satisfy (after A- homogenization as in Lemma IXTj) 
different systems in the family Hj, i = 0, . . . , m^~^ — 1 according to (/). Each system 
is satisfied precisely by m of these homogenized functions, which we denote by zi{a). 
The functions zi{a) also satisfy the corresponding algebraic equation: 

(4.8) amz"" + e'^aiz""' + e''a2z"'' + ...+ e'^anz""" + a^+i = 0. 

In particular, when I = (ii, . . . , i„) with ik = jmk mod (m) for all A; = 1, . . . , n and 
some fixed j = 0, . . . , m— 1, the function zj satisfies the differential yl-hypergeometric 
system recalled in Section 01 and the algebraic equation (jLip . Recall that yi{x) = 
zi{l,Xi, . . . ,Xn, —1) as in Theorem 13.31 Changing a little bit the notation of the 
homogeneous roots as in (|4.4p . in order to prove that the function Xc is a solution 
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of the Mellin system, it is enough to prove that the homogenized function 

(4.9) Xc ■■= E ^^'^(^!'^ + ■ ■ ■ + log(4'^)) 

fc=l 

satisfies the lattice hypergeometric system Tig. Then, it is possible to prove that 
each of the functions 

(4.10) a^^^ := zf^ log(zf ^) + ■ • • + log(4'^) 

is annihilated by the toric operators in the corresponding torus translate of the A- 
hypergeometric system as in the proof of Proposition 3.2 in 4J. Moreover, the fact 
that c lies in R allows one to prove that the whole sum Xc satisfies the homogeneity 
equations. 

The theorem will be proved if we show the following statement: If a linear com- 
bination x\ '■= J2i=i '^iXciii is algebraic, where Aj G C and c[l], . . . ,c[s] e R, then 
J2i=i = 0. Note that this implies in particular (since the function is alge- 
braic) that the dimension of the linear span of the functions {xc '■ c E R} equals 
the dimension of R, from which we deduce not only that the sum 5 © F is direct 
but also that it equals the full solution space of the Mellin system. 

In order to prove the above statement, we observe that if x\ is an algebraic 
function of (xi, . . . , Xn) then, the homogenized function 

i=l k=l 

is an algebraic function of a = (om, ami, • • • 5 ^mn, dn+i)- But now, as in [U Th. 3.5], 
we fix a point a with ja^+il sufficiently small relative to the absolute values of the 
other coordinates and consider the loop 7 defined by 

Analytic continuation around 7 will return all the local roots in a suitable neighbor- 
hood of a to their original values, while there is a subset S^''^ containing m„ < m 
of the roots, such that the logarithms of the roots in S^''^ will be increased by 27rz. 
Then 7*(a(^-))-«('^) equals 2TTi times the sum of the m„ roots in S^''\ which cannot 
be zero since R = R' . We denote this sum of roots by Zg^k). Therefore 

i=l k=l / i=l k=l 

which, after dividing by 27ri, equals 

Ea. E = E E^^^w^'^ 

i=l fc=l fc=l \i=l / 

Since z^g'^ is not the complete sum of all the roots corresponding to one equation, 
the only way in which a finite iterate of the action of 7 returns X\ to its original 
value is when Yl!i=i Ajc[z]*-'^^ = for all fc, i.e., Yll^i \c[i\ = 0, as asserted. □ 
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We also deduce: 

Corollary 4.5. The monodromy representation of the Mellin system M.'Jji is always 
reducible. 

Proof. Suppose that n> 1. By Theorem 12.11 the dimension of the space of analytic 
solutions to ()1.3|1 equals m". Since the solutions to the algebraic equation ()1.2|) span 
a linear space of dimension at most m, there is a nontrivial invariant subspace in 
the solution space of (|1.3|) . If n = 1 and mi < m — 1, then the solutions of the 
algebraic equation form a nontrivial invariant subspace of dimension m — 1. Finally, 
for n = 1 and nti = m — 1, the one-dimensional linear space spanned by the rational 
solution Xi is invariant under the action of the monodromy. □ 



5. The case of dimension one 

In this section we further investigate the special case of the Mellin system of 
equations associated with the trinomial algebraic equation 

(5.1) I/*" + xy'"! - 1 = 0. 

In this case the Mellin system ()1.3|) consists of a single ordinary differential equation 

of order m, namely the equation 

(5.2) 

''mi — 1 m— mi — 1 \ 

JJ(mi^ + mfc + l) JJ ((m-mi)^ + mfc- 1) - (-l)'"im'"-^ U(x) = 0, 

fc=0 A:=0 / 

where 6 = x-^. Let m = da, mi = db, where d = GCD(m, mi). The discriminant 
of the algebraic equation (j5.1|l (computed with respect to y) is given by x"'b^{a — 
bY~^ — (— l)^a". For d = 1, this expression coincides with the leading coefficient in 
the differential equation ()5.2j) . Throughout this section, we will call the differential 
operator in the left-hand side of ()5.2p the Mellin operator and denote it by M{m, mi). 
Special instances of this operator are considered in Examples 12 . 315 . 1 15 . 31 and 14.11 

Example 5.1. Consider the cubic equation 

(5.3) y^ + xy^ -1= 0. 

The corresponding Mellin equation can be written in the form 

(5.4) (4x^ - 27)y"' + ISx^y" + Axy' - 4?/ = 0. 

A basis in the solution space of the differential equation (|5.4|) is given by the functions 



yi{x)=x, y2{x) = Y 108 - Sx^ + 12V81 - 12x3, y^{x) = x^ / y2{x) , 

or by the roots of the algebraic equation ()5.3p . 

Notice that after multiplication with x^ the differential operator defining the equa- 
tion ()5.4|) can be factorized as follows: 

x^ f (4x^ - 27)-^ + ISx^-^ + 4x-^ - 4 
V dx'^ dx^ dx 
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( {4x^ - 27x)^ + {Ux' + 27)4- + 4x'l ( x^- - l] ■ 
\ dx^ ax J \ dx J 

This factorization makes the rational solution yi{x) = x obvious. 

More in general, assume mi = m — 1. Since the roots yi{x), . . . ,ym{x) of the 
trinomial satisfy the equation ()5.2j) and sum up to x, it follows that the operator 6 — 1 
(which is the annihilator of x) can be factored out of M{m,m — 1). This factorization 
is described in the following proposition. 

Proposition 5.2. For mi = m — 1, the only nontrivial decomposition of the Mellin 
operator defining the equation i5.^} in the Weyl algebra of linear differential opera- 
tors with polynomial coefficients is given by 

(m—2 m—1 \ 

x"" Yl ((m - 1)^ + mk + l) + (-m)*^ 9 Y[{9 - k)] {9 - 1). 
k=0 k=2 / 



The proof of this proposition is analogous to the proof of Proposition 15.41 below. 
Assume now that m > mi + 1 and d = GCD (m, mi) = 1. 

Example 5.3. Consider the cubic equation 

(5.6) y^ + xy -1 = 0. 

The corresponding Mellin equation can be written in the form 

(5.7) {Ax^ + 27)y"' + ISx'^y" + lOxy' - 2y = 0. 

Since the coefficient by y"^ in ()5.6p is zero, the solutions to ()5.6|) generate a two- 
dimensional vector space and thus there exists a solution x to the differential equa- 
tion (|5.7|) which is not a solution to the algebraic equation (j5.6j) . The following 
functions give a basis of solutions to the Mellin equation ()5.7j) : 

zi{x) = \/l08 + 12 Vl2x3 + 81, Z2{x) = x/zi{x), 

Xix) = yi log(yi) + y2 log(?/2) + Z/s log(?/3), 

where yi{x) , y2{x) , y^^x) are the three (linearly dependent) solutions to ()5.(i|l . 

Notice that the differential operator defining the equation ()5.7|) can be factorized 
as follows: 

(4x3 + 27)-^ + 18x^4^ + lOx-^ - 2 = f (4^^ + 27)-^ + Gx^-^ - 2x 
dx'^ dx^ dx dx \ ax^ dx 

The functions zi, Z2 give a fundamental system of solutions to the second-order factor 
in this decomposition. This operator is described in Pj. 

It turns out, that the Mellin operator corresponding to mi = 1 allows a similar 
decomposition. 
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Proposition 5.4. For mi = 1, the only nontrivial decomposition of the Mellin op- 
erator defining the equation i5.^) in the Weyl algebra of linear differential operators 
with polynomial coefficients is given by 

d ( \ 

(5.8) M(m, 1) = — X TT ((^ - 1)^ + "^'^ - 1) + 7 • 

dx \ dx^'^ ] 

\ fc=0 / 

Proof. The validity of ()5.8|) follows from the Weyl algebra identity 9 + 1 = ^ o 
X. Here 'o' denotes the composition of differential operators in the Weyl algebra. 
Let us denote the second factor in the right-hand side of ()5.8|) by Maig. We would 
like to show that the differential equation M.^\^y = has irreducible monodromy. 
Multiplying this equation with a;™~^, using the identity (j2.7p and making the change 
of variables t = x™, we arrive at the equation 

(5.9) ((,„ - 1)"- 1 n [e. + + n {o, -i))y- 0. 

where 9t = t-^. By Proposition 3.3 in [2j the monodromy of the differential equa- 
tion (|5.9j) is reducible if and only if {mi — l)/m{m — l)+j/m& Z, for some i,j = 
0, . . . , m — 2. This condition can only be satisfied when {mi — l)/m{m — l) +j/m = 1. 
Clearing the denominators and considering the obtained relation modulo m we con- 
clude that it can never be satisfied. Thus the monodromy of ()5.9|) is irreducible and 
hence Maig does not admit any further decomposition. □ 

In case d > 1 the Mellin operator always factors. 

Proposition 5.5. There exist d linear operators Qo, . . . , Qd-i of order ^ with co- 
efficients in C{x) such that M{m, mi) = Qq . . . Qd-i- 

Proof. According to Corollary 12.51 and the discussion in Section HI all solutions to 
M{m,mi) are algebraic (so its associated Galois group coincides with the mon- 
odromy group) and a basis of solutions is given by the functions 

{/,, ,(x) := e^yprie'^^^+'x), j = 0, . . . ,m/d - 1, k = 0, . . . ,d - 1}, 

where e is a primitive m-root of 1 and ypr is the principal solution described in ()2.1|) . 
For each k = 0, . . . ,d — 1, the m/d functions {fj,k{x), j = 0, . . . , ma-i} are the 
(linearly independent) roots of the algebraic equation 

y'^ + e'^xy"'' -1 = 0, k = 0,...,d-l. 

Then, they span a monodromy invariant subspace Vk of the solutions V to M{m, mi) 
of dimension m/d, which gives a right factor of order m/d of the Mellin operator 
whose solution space equals V^. Moreover, we have the direct sum decomposition 
V = Vo©- ■ -©Vrf-i. Then, it follows by the results in [T3] that there exist operators 
Qo, ■ ■ ■ ,Qd-i (with Qi equivalent to Li in the sense of [TS|) such that M{m,mi) 
decomposes as M{m, mi) = Qq . . . Qd-i- D 
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Example 5.6. Let n = 1, m = 4, mi = 2. The Mellin equation associated with the 
algebraic equation + xy"^ — 1 = is given by 

(5.10) (16x^ - 256)?/(^) + WOx'^y'" + 360x^y" + I20xy' - I5y = 0. 

In accordance with Corollary 12 .51 a basis in the solution space of ()5.1()j) is given by 
solutions to the algebraic equations y^ + xy"^ — 1 = and y^ + ixy"^ — 1 = 0. After 
a suitable normalization this basis is given by the functions \/ —x ± Va;^ ± 4. The 
solution space of ()5.10|) splits into two 2-dimensional invariant subspaces spanned 

by a/— a; ± Va;^ + 4 and \/ —x ± a/z^ — 4. This decomposition of the solution space 
corresponds to the factorization of the differential operator defining ()5.10p : 

(74 ^3 £. ^ 

(16x^ - 256)— + leOx^— + 360x2— ^ i20x— - 15 = 

xajJu KAjtXj \AjtXj LtitXj 

, o ^ d? d \ f , n X c?2 d 

4x2 - 16 — + 20x— + 15 4x2 + 16 — + 4x— - 
dx^ dx I \ dx'^ dx 



Example 5.7. Let n = 1, m = 6, mi = 2 and consider the algebraic equation 
(5.11) + x?/2 - 1 = 0. 

The Mellin equation associated with 1)5.111) is given by 
(1024x^ - 46656)?/('') + 27648x^?/(^) + 242816x^?/(^) + 



(5. 12) 818944x^2/'" + 955780x2?/" ^ 236180x?/' - 6545y = 0. 

A basis in the space of solutions to ()5.12p is given by {/(e^x), i = 0, . . . , 5}, where 
e = e^"^/^ and /(x) is the principal solution to 1)5.11)) . The differential operator 
defining the equation ()5.12)) admits the following factorization: 

d^ d^ d^ 

(1024x^ - 46656)— + 27648x^— + 242816x^— + 

^ ' dx^ dx^ dx^ 



818944x^— + 955780x2— ^ 236180x— - 6545 
dx-^ dx^ dx 



d^ d^ d 

(32x^ - 216)— + 432x2— + 1526x— + 1309 

dx'^ dx^ dx 



d^ ^2 d 

(32x^ + 216)— + 144x2— + 86x— 
ctx'^ dx^ dx 
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